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Abstract. In most experiments on the Casimir force the comparison between 
measurement data and theory was done using the concept of the root-mean-square 
deviation, a procedure that has been criticizcd in literature. Here we propose a special 
statistical analysis which should be performed separately for the experimental data and 
for the results of the theoretical computations. In so doing, the random, systematic, 
and total experimental errors are found as functions of separation, taking into account 
the distribution laws for each error at 95% confidence. Independently, all theoretical 
errors are combined to obtain the total theoretical error at the same confidence. Finally, 
the confidence interval for the differences between theoretical and experimental vàlues 
is obtained as a function of separation. This rigorous approach is applied to two recent 
experiments on the Casimir effect. 



PACS numbers: 12.20.Fv, 12.20.Ds, 07.05.Kf 



1. Introduction 

Today the Casimir effect is being actively investigated not only theoretically but also 
experimentally. Historically the first measurement of the Casimir force between metals 
was performed in 1958 pQ and confirmed the existence of the force with an uncertainty 
of about 100%. In the following dècades the experimental output was painfully low and 
only one experiment with metal test bodies was made [2j (see Ref. [3] for a review). In 
the last few years many measurements of the Casimir force have been performed using 
torsion pendulums, atòmic force microscopes, micromechanical torsional oscillators, and 
other laboratory techniques [4-16]. Most authors (see Refs. [1-14]) have used the 
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concept of the root-mean-square deviation between experiment and theory to quantify 
the precision of the measurements. However, for strongly nonlinear quantities, such as 
the Casimir force which changes rapidly with separation distance, this method is not 
appropriate because it may lead to different results when applied in different ranges of 
separations. This was emphasized in Ref. jH] although no better method was suggested. 

The present paper contains the comparison analysis of the precision and accuracy in 
two recent experiments ITü] using rigorous methods of mathematical statistics. The 
distinctive feature of our approach is that both total experimental and total theoretical 
errors are determined independently of one another at some accepted confidence level. 
Then, the absolute error of differences between calculated and measured vàlues of the 
physical quantity is found at the same confidence as a function of separation, serving as 
a measure of the precision in the comparison of experiment and theory. 

2. Determination of the experimental errors 

2.1. Random errors 

In experiment the Casimir pressure between two Au coated parallel plates was 
determined dynamically by means of a microelectromechanical torsional oscillator within 
the separation region from 160 to 750 nm. In experiment the Casimir force was 
measured between a Si plate and a large Au coated sphere using an atòmic force 
microscope within the separations from 62.33 to 600.04 nm. In our error analysis we use 
the notation II (z) which denotes either the measured Casimir pressure P exp (z) or force 
F exp (z) as a function of separation z between the test bodies. 

Usually several sets of measurements, say n, are taken within one separation region 
(z a , Zb). This is done in order to decrease the random error and to narrow the confidence 
interval. In Ref. ^2] n = 14, and in Ref. Jïï] n — 65. Each set consists of pairs [z%, n(zj)] 
where 1 < i < w = (288 -l· 293) in Ref. [IHJ and 1 < i < w = 3164 in Ref. [16] . All 
measurement data can be represented by pairs [zij, 11(2^)] where 1 < j < n. Generally 
speaking, separations with fixed i but different j may be different (this was the case in 
Ref. [IS])- For such measurement results it is reasonable to divide the entire separation 
range (z a , zj,) into partial sub intervals of length 2Az, where Az is the absolute error in 
the measurement of separations equal to 0.6 nm and 0.8 nm in Refs. |15[ll6j. respectively. 
In so doing, each subinterval k contains a group of m k points Zy = z kh 1 < l < m k 
(in Ref. mk ranges from 3 to 13). Inside each subinterval all points z^i can be 
considered as equivalent, because within the interval of width 2Az the value of absolute 
separation is distributed uniformly. The mean and the variance of the mean of the 
physical quantity II for the subinterval k are defined as 



If all z^ = Zi, i.e., the same in different sets of measurement (as in Ref. [HI|), the 




(1) 
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mean and the variance of the mean at each point Zi are obtained more simply 

(2) 



i n i n 



Direct calculation shows that the mean vàlues fïfc, ÍT are uniform, i.e., change 
smoothly with the change of k, i. The variances of the mean, Sfi k , Sfj., are, however, not 
uniform. To smooth them, we have used a special procedure developed in mathematical 
statistics ^7J QB]. At each separation z , in order to find the uniform variance of a 
mean, we consider not only one subinterval containing zq but also several neighboring 
subintervals from both sides of z (4 or 5 in Ref. [15J) or about 30 neighboring points in 
Ref. [TSj. The number of neighboring subintervals or points is denoted by N. Then the 
smoothed variance of the mean at a point z is given by [T7[ HH] 

N 



s ü( z o) — max 



n E x l s k 



k=l 



(3) 



where are the statistical weights. The maximum in Eq. J3J) is taken over two sets of 
coemcients, = 1/N, and \ k = l/(c^ J2fLi cï l ) where the constants c$ are determined 
from Sg : s\ 2 : . . . : s\ — C\ : Ci : . . . : Cat. Note that max in Eq. (JHJ) leads to the 
most conservative result, i.e., overestimates the random error. Finally, the confidence 
interval at a confidence probability f3 takes the form 

n(2b) - A rand n(z ), n(^ ) + A rand n(^ )] , (4) 

where the random absolute error in the measurement of the quantity II at a separation 
z is given by 

A rand n(,2o) = s n (z )t (l+m (mmm k - 1). (5) 

Here the value of t p (f) can be found in tables for the Student's í-distribution. For 
example, in the experiment |16j. minmfc = n = 65. Thus, for (3 = 0.95, we have 
t p (f) = 2.00 and A rand F cxp = 3.0 pN is independent of z . 

The computational results for the relative random errors <5 rand n = A rand n/|n| in 
the experiments f^lE] at 95% confidence are shown in columns labeled (a) in Table 1, 
as the functions of separation. As is seen from column two in Table 1, in the experiment 
[TH] the relative random error of the Casimir pressure measurements is equal to 1.5% at 
z = 160 nm, then it quickly decreases to 0.4% at z — 350 nm, and then increases with 
further increase of separation. This is explained by the fact that the absolute random 
error in Eq. © takes a maximum value at the shortest separation and monotonically 
decreases with the increase of separation until z = 400 nm [TS] . At larger separations 
A rand P exp is practically constant and the increase of 5 rand P exp is explained by solely 
in terms of the decrease of the Casimir pressure magnitude. In the experiment Jü] 
(column 8 in Table 1) the absolute random error is only 0.78% at the shortest separation 
z = 62.33 nm, and quickly increases with separation due to the decrease of the Casimir 
force. 
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Table 1. Relative errors (%) in experiments ^3 Ej: random errors 

f and F cx P, ,f and^cxp ( a ). sys tematic errors S syst P cxp , S syst F ex P (b); total experimental 
errors S tot P CK P, S tot F cx P (c); theoretical errors S a P th , S a F th (d); total theoretical errors 
S tot P th : S tot F th (e). Columns labeled (f) contain S p /|P cx p| and E F /\F e * p \ (see text). 



Experiment of Ref. |15| Experiment of Ref. |16j 



z (nm) 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


62.33 














0.78 


0.31 


0.87 


0.55 


3.5 


4.0 


70 














1.1 


0.42 


1.2 


0.56 


3.2 


3.7 


80 














1.6 


0.60 


1.7 


0.56 


2.8 


3.7 


90 














2.1 


0.84 


2.4 


0.56 


2.6 


3.9 


100 














2.9 


1.1 


3.2 


0.56 


2.4 


4.4 


120 














4.7 


1.8 


5.3 


0.56 


2.0 


6.2 


140 














7.3 


2.8 


8.1 


0.57 


1.8 


9.1 


160 


1.4 


0.15 


1.4 


0.56 


1.6 


2.4 


10 


4.1 


12 


0.58 


1.6 


13 


170 


0.59 


0.15 


0.59 


0.56 


1.6 


1.9 


12 


4.9 


14 


0.58 


1.6 


15 


180 


0.57 


0.15 


0.57 


0.57 


1.5 


1.8 


15 


5.7 


16 


0.58 


1.5 


18 


200 


0.55 


0.16 


0.56 


0.57 


1.4 


1.7 


20 


7.7 


22 


0.59 


1.4 


23 


250 


0.48 


0.20 


0.54 


0.58 


1.2 


1.5 


37 


14 


41 


0.61 


1.3 


42 


300 


0.44 


0.31 


0.59 


0.59 


1.1 


1.4 


62 


24 


69 


0.64 


1.2 


70 


350 


0.40 


0.50 


0.72 


0.61 


1.0 


1.4 


96 


37 


107 


0.67 


1.1 


108 


400 


0.56 


0.80 


1.1 


0.62 


0.98 


1.6 














500 


1.3 


1.80 


2.5 


0.66 


0.91 


2.9 














600 


2.9 


3.80 


5.4 


0.70 


0.88 


5.4 















2.2. Systematic errors 

In each of the experiments [To'l lTïï] there are several absolute systematic errors òtj st lí{z) 
and respective relative systematic errors 5^ yat H(z) = A h / st II(2:)/|ri(,2:)| where 1 < i < J. 
Systematic errors are the random quantities characterized by a uniform distribution. 
Because of this, the total systematic error is JÜ] 




5^Tl{z) = min £ ^(z), £ ^(z) , (6) 



where (3 is the confidence probability, and k^p is a tabulated coefhcient |T5] . The same 
rule is also vàlid for the absolute systematic errors. 

In the experiment ^Hj there are J = 2 main systematic errors: 

^systpexp = 5R = ^ ^ St p-P( 2 ) = S(üü r - U ) = ^ (7) 
R \üü r — üüq\ 

where R = (148. 7± 0.2) fim is the sphere radius, u r and u> are the resonant and natural 
angular freqüències of the oscillator, respectively (the former is separation-dependent). 
u)o = 2n x 702.92 Hz was determined so precisely that its error does not contribute to 
the results, and the error of the resonant frequeney is Au r = 2ir x 6mHz. Using the 
value /cq.95 = 110 an d utilizing Eq. (jüj) one obtains the total systematic errors given in 
column 3 [labeled (b)] in Table 1. 
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The experiment [16J contains the following J = 4 systematic errors: AY st -F exp ~ 
0.82 pN due to force calibration; A s 2 yst F exp w 0.55 pN due to noise when the calibration 
voltage is applied to the cantilever; A S J S F exp ~ 0.31 pN due to the instrumental 
sensitivity; and A s / st .F exp ~ 0.12 pN due to the restrictions on computer resolution 
of data. Combining these errors using the analog of Eq. © with fcg 4 g 5 — 112, we obtain 
A syst F c xp = i.i7 p N. The respective relative errors ó syst F cxp = A syst F cxp /|F exp | are 
shown in column 9 in Table 1. Comparing columns labeled (b) in Table 1, we conclude 
that in both experiments the relative systematic error increases as the separation 
increases. The magnitudes of the systematic errors are smaller in the experiment of 
Ref. US]. 



2.3. Total experimental error 

To find the total experimental error in the measurements of 11(2), one should combine 
the random and systematic errors obtained above which are described by a normal 
(or Student) distribution and a combination of uniform distributions, respectively. To 
be very conservative, we assume that the systematic error is described by a uniform 
distribution (other assumptions lead to smaller total error). Different methods for 
combining random and systematic errors are described in Ref.|19j. Here we use one 
based on the value of the quantity r(z) = A syst n(z) / Sjj(z) . According to this method, 
at all z where r(z) < 0.8 the contribution from the systematic error is negligible and 
A tot n(z) = A rand n(z) at 95% confidence. If r(z) > 8 is vàlid, the random error is 
negligible and at 95% confidence A tot n(z) = A syst IT(z). In the separation region where 
0.8 < r(z) < 8, the combination of errors is performed using the rule 

A tot n(^) = q p (r) [A rand n(^) + A syst IT(z)] , (8) 

where the coefficient qp{r) with (3 = 0.95 varies between 0.71 and 0.81. Being 
conservative, here we use qp{r) = 0.8 in all calculations. 

Table 1 [columns 4 and 10 labeled (c)] contain the total experimental error of the 
Casimir pressure and force measurements in the experiments |15| ITïïj , respectively. As 
seen in column 4 of Table 1, in the experiment ^H] at z = 160 nm the total experimental 
error is equal to 1.4%, but in a wide separation range from 170 to 300 nm, it is practically 
flat and within the range from 0.54 to 0.59%. Even at z — 600 nm it is equal to only 
5.4%. In the experiment JH)| (column 10 in Table 1) the smallest total experimental 
error of 0.87% is achieved z = 62.33 nm and increases up to 5.3% at z = 120 nm. This 
is mainly due to the large contribution of the random errors. 



3. Determination of the theoretical errors 



The theoretical vàlues of U.(z) (both the pressure and force) are computed using the 
Lifshitz formula (see, e.g., Ref. [3J which takes into account the effects of finite 
conductivity and nonzero temperature. The Lifshitz formula contains the reflection 
coefficients at imaginary Matsubara freqüències. At zero Matsubara frequency these 
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coeficients are expressed in terms of the Drude dielectric function (the Drude model 
approach |2"Ü1 I2~T| ) or in terms of the Leontovich surface impedance (the impedance 
approach [221 ES])- At nonzero Matsubara freqüències both approaches use the tabulated 
optical data extrapolated to low freqüències by the imaginary part of the Drude dielectric 
function. In Refs. [2U EE] the reflection coefficients at all Matsubara freqüències were 
expressed using the free electron plasma model (the plasma model approach). 

One error in the theoretical computation arises from sample to sample variations 
of the optical data for the complex index of refraction. Usually these data are not 
measured in each individual experiment, but are taken from Tables. In Ref. [2] it was 
shown that variation of the optical data for typical samples leads to a relative theoretical 
error 5\II til (z) in the computed Casimir pressure or force that is no larger than 0.5%. 
Being conservative, we set SiU th (z) = 0.5% at all separations. Strictly speaking, there 
may occur rare samples with up to 2% deviations in the Casimir pressure or force at 
short separations. If this happens, the theoretical vàlues come into conflict with the 
experimental data. Such deviations must be considered not as an error (they can only 
diminish the magnitudes of the pressure or force) but as a correction. The validity of the 
hypothesis on the presence of such types of corrections can be easily verified statistically. 

Another theoretical error is caused by the use of the proximity force theorem |2b'| . 
(This is the name given by the authors of Ref. [2H]; some other authors, e.g. in Ref. |27j . 
prefer to use the name "proximity force approximation" to underline the approximate 
character of the equality proposed in Ref. [2E]-) In the experiment |T5] it is applied to 
express the effective Casimir pressure between two parallel plates through the derivative 
of the force acting between a sphere and a plate. In the experiment Jïï] the bàsic result 
for the force is obtained using the proximity force theorem. The upper limit of error 
introduced by this is ó 2 Tl th (z) = z/R [3] (see also Refs. [23121] where the same estimation 
was confirmed for the case of a massless scalar field). 

Both errors 8jjl are described by a uniform distribution and in this sense can be 
likened to systematic errors. They are combined by using Eq. (JïïJ) with J = 2 leading 
to the vàlues 5oIl th presented in columns 5 and 11 in Table 1 [labeled (d)] for the 
experiments fSlQBl, respectively. As is seen from these columns, the errors 5oIl th (z) 
depend only slightly on separation and take similar vàlues between 0.55 and 0.70%. 

In addition to the major theoretical errors 5iU th (z), there exist other uncertainties 
in calculations which are not taken into account in the Lifshitz formula. Some of them 
were shown to be negligibly small (like the contributions from patch potentials, nonlocal 
effects and finite sizes of the plates |141 115j ). As to the contribution from the surface 
roughness, it was calculated using the atòmic force microscope images of the interacting 
surfaces and taken into account as a correction [14-16]. This is why these factors do 
not contribute to the balance of theoretical errors. 

There is one more error which can be considered together with the theoretical errors 
if one is going to compare the experimental and theoretical vàlues of Tí(z) ^3^21- This 
arises from the fact that z is determined experimentally with an error Az (see Sec. 2.1), 
and this error results in the additional uncertainties S^Il (z) in computations. Bearing 
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in mind the leading theoretical dependences of the pressure and force on separation, 
we obtain 5 3 P th (z) = AAz/z in Ref. [T5] and 5 3 F th (z) = AR/R + 3Az/z in Ref. 
Taking into account that the combined random quantity SqH (z) may be distributed 
nonuniformly, we combine it with 5 3 U th (z) using Eq. (JSJ) and obtain the total theoretical 
error ó tot Tl th (z) at 95% confidence. The vàlues of 5 tot n th (z) are presented in columns 
6 and 12 in Table 1 [labeled (e)] for the experiments respectively. For both 

experiments they monotonically decrease with separation and take the largest vàlues at 
the shortest separation. The significant increase of the total theoretical error in columns 
labeled (e) compared to those labeled (d) is due to the additional error 5 3 Il th (z). 

4. Comparison between experiment and theory 

4-1. Measure of agreement between experiment and theory 

In Secs. 2.3 and 3.2 we have obtained the total experimental and theoretical errors 
at 95% confidence for both the Casimir pressure and force. Now we consider the new 
random quantity P th (z) — P exp (z) [or F th (z) — F exp (z)] and determine the absolute error 
of this quantity, Ep tF (z), at 95% confidence using the composition rule (jüj) with J = 2 

E P (z) = min |A tot P th (^) + A tot P cxp (^), 1.1 ^[A tot P th (z)] 2 + [A ot P cxp (z)] 2 ^ (9) 

(the same equation is vàlid for the force). Note that in Eq. (JÜJ) the conservative value 
of k^l 5 = 1.1 is used as for two uniform distributions (otherwise it would be smaller). 

The confidence interval for the quantity P th (z) — P cxp (z) at 95% confidence is 
given by [—Ep(z), E, P (z)) and the mean vàlues (P th (z) — P cxp (z)) or (F th (z) — F exp (z)) 
must belong to this interval or its analog for the force with a 95% probability. The 
vàlues of 'Ep(z)/\P exp \ and 'Ep(z)/\F exp \ are given in columns 7 and 13 in Table 1 
[labeled (f)] for the experiments JEIEI], respectively. They characterize the sensitivity 
of the experiments fü[ to the differences between theory and experiment at 95% 
confidence. For example, in Ref. ^H] theory is in agreement with experiment at a 
separation z = 400 nm if \P th (z) — P cxp (z)\ does not exceed 1.6% of \P exp (z)\. 

4-2. Comparison of experiment and theory in the determination of the Casimir 
pressure between Au plates in Ref. 115^ 

Experiment ^3] is rather sensitive and can be compared with different theoretical 
approaches to the calculation of the Casimir pressure. The main results are presented in 
Table 2 where the second and third columns labeled (a), (b) contain the half-width Hp(z) 
of the confidence interval at 95 and 99% confidence, respectively. In columns 4-7 labeled 
(c), (d), (e) and (f) the results for the mean differences (P th (z) — P cxp (z)) are computed 
using the impedance [13 122] and the plasma model [13 12H I2ü] approach at T = 300 K, 
the optical data in the Lifshitz formula at T = [Hj, and the Drude model approach 
at T = 300 K [2H1I2I1, respectively. To avoid confusion, recali that in column (c) the 
zero-frequency contribution to the Lifshitz formula is computed using the Leontovich 



Comparison between experiment and theory in Casimir force measurements 



8 



Table 2. Comparison of the experiments of Refs. |15l I16| with theory. Columns (a) 
contain the absolute errors Ep t F(z) of the pressure (mPa) and force (pN) differences 
at 95% confidcncc. Column (b) contains the same quantity for the pressure at 99% 
confidence. Other columns contain the mean vàlues (P th (z)— P exp (z)) in mPa [and also 
(F th (z) — F exp (z)) in pN for the last column labeled (e)] computed using four different 
approaches: the impedance (c) and the plasma model (d) approach at T = 300 K; 
the optical data in the Lifshitz formula at T — (e); the Drude model approach at 
T = 300 K (f). 

Experiment of Ref. |15| Experiment of Ref. |16) 



z (nm) 


fal 


fbl 


fel 


fdl 


fel 


ffl 


fal 


fel 


62.33 














15.2 


-0.5 


70 














10.4 


3.0 


80 














7.1 


3.6 


90 














5.4 


1.0 


100 














4.5 


2.0 


120 














3.9 


-0.15 


140 














3.8 


0.02 


170 


17.2 


39.8 


2.01 


13.0 


3.87 


18.8 


3.7 


-0.82 


180 


13.4 


31.0 


-0.74 


7.54 


1.24 


14.4 


3.7 


-0.48 


200 


8.59 


19.8 


-1.21 


5.3 


0.63 


11.0 


3.7 


-0.31 


250 


3.34 


7.72 


-0.31 


1.3 


0.93 


7.09 


3.7 


-0.84 


300 


1.59 


3.67 


0.34 


0.6 


1.12 


5.07 


3.7 


0.46 


350 


0.89 


2.06 


0.38 


0.39 


0.80 


3.58 


3.7 


0.27 


400 


0.63 


1.46 


0.28 


0.20 


0.68 


2.59 






500 


0.49 


1.13 


0.11 


0.05 


0.32 


1.37 






600 


0.46 


1.06 


0.08 


0.04 


0.17 


0.82 






700 


0.46 


1.06 


0.02 


-0.01 


0.08 


0.51 







impedance in the region of infrared òptics. At all other Matsubara freqüències the 
impedance is obtained using the tabulated optical data. Comparing columns 4-6 and 
columns 2,3, we conclude that the impedance approach, the plasma model approach and 
the Lifshitz formula at T = are consistent with the measurement data. At the same 
time, by comparing columns 2,3 with column 7 we find that the Drude model approach 
is excluded by experiment at 95% confidence within the separation range from 170 to 
700 nm, and at 99% confidence from 300 to 500 nm. The physical reasons for the failure 
of the Drude model approach and the advantages of the Leontovich impedance are 
discussed in Refs. |221 Efl E2| • 

4-3. Comparison of experiment and theory in measuring the Casimir force between an 
Au sphere and a Si plate in Ref. fTpy 

Experiment [Ej is the first demonstration of the Casimir force between a metal and 
a semiconductor performed at shorter separations than in experiment jT^j. For this 
reason it cannot be used to discriminate among different theories. In column 8 in Table 
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2 labeled (a) the vàlues of E F (z) for the force at 95% confidence are given. Column 9 in 
Table 2 labeled (e) contains the vàlues of (F th (z) —F oxp (z)) computed using the Lifshitz 
formula at T = and tabulated optical data for Au and Si. The comparison of these 
columns shows that the theory at T = is in a very good agreement with experiment. 

5. Conclusions 

From the above, several conclusions can be reached: 

- A new method for data processing and comparing theory with experiment for the 
Casimir effect has been presented based on rigorous results of mathematical statistics 
with no recourse to the previously used root-mean-square deviation; 

- The distinguishing feature of this method is the independent determination of 
the total experimental and theoretical errors and of the confidence interval for differences 
between calculated and measured vàlues at a chosen confidence probability; 

- The developed method is conservative and guaranties against underestimation 
of errors and uncertainties. It was applied to two recent experiments measuring the 
Casimir pressure and force in different configurat ions; 

- We have demonstrated that the approaches based on the vanishing contribution 
of the transverse elèctric mode at zero frequency (e.g., the Drude model approach) are 
excluded by experiment at 99% confidence, whereas the three traditional approaches to 
the thermal Casimir force are consistent with experiment. 
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